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This paper presents the constraints on the time variation of the fine structure con- 
stant at recombination relative to its present value, A a = (a rec — ao)/ Q (b ob- 
tained from the analysis of the WMAP-3yr Cosmic Microwave Background (CMB) 
data, with an additional prior on the Hubble expansion rate from HST Hubble 
Key Project. I found out that -0.039 < A a < 0.010 at 95% C.L., which brings a 
30% improvement to the previous limits from WMAP-lyr data. The correspond- 
ing recombination redshift, z rec = 1075lgg, shows a delayed recombination epoch 
compared with the results from WMAP-lyr data. 
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1 INTRODUCTION 



With the rele ase of WMAP-3yr temperature (IHinshaw et al. L 120061 ) and polar- 



ization data (jPage et al. 1 . 120061 ) . the standard ACDM cosmolo gical model has 



passe d another rigorous set of tests, being validated once again (ISpergel et al. 



20061 ) . In this era of precision cosmology, the early Universe is playing an in- 
creasingly more important role as a laboratory where physics not accessible 
on Earth can be tested and verified. Thus, cosmology became an excellent 
tool to answer to many questions of fundamental physics, among which the 
time variation of fundamental constants of nature, predicted a s observational 



consequences of some currently preferred unification theories ( Isee e.g. Uzan 



20031 : Martins I . 12002 ). 
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One of the most exciting case is the fine structure constant a, possible time 
variations being already reported from geophysical and astronomical observa- 
tions. 



From the Oklo natural nuclear reactor in Gabon, the relative time variation 



of a has been obtained as —0.9 x 10 7 < A„ < 1.2 x 10 7 at z 



0.1 



( IDamour fc Dyson I . Il996l ) 



Terrestrial laboratory measurements give the constraint Aa/a < 1.4 x 10 14 
based on comparisons of rates b etween clocks with dif ferent atomic numbers 



during 140 days of observations (IPrestage et al. 1 . 119951 ). 



While the results presented above are only upper limits of A a , the estima- 
tions of its actual value have been firstly computed from spectral analysis 
of high-redshift quasar absorption systems. Some of the reported values are 
A„ = (—4.6 ± 5.7) x 10~ 5 for a redshif t range z ~ 2 — 4 (IVarshalovich et al 



l200lh . A a = (-1.09 ± 0.36) x 10~ 5 f lWebb et all l200lh . these determina- 



tions being a confirmation of the expectation that the fine structure con- 
stant is a nondecreasing function of time. From this point of view, the re- 
su lt (a(z 9 ) - a(z- i) )/a n = (5.43 ± 2.52) ppm for z x = 1.15 and z 2 = 1.84 



( ILevshakov et al. 1 . 120061 ) . is controversial. In fact, there are also many other 



contradictory results which, perhaps, may be validated if the fine structure 
constant would not only be time dependent but also had a spatial variation. 

The CMB presents the advantage of probing directly the decoupling epoch 
[z ~ 1 100) when the a variations from its current value are expected to be 
more important and therefore, more easy to test. A value of a different from ao 
at the epoch of recombination changes the ionization evolution at that epoch, 
inducing effects on the CMB anisotropics observed today. 

In the last years several groups have analyzed the variation of fine structure 
constant at recombination using CMB measurements. A lower value of a at 
recombination epoch compared with its present value, was taken into account 
as a possibility to solve the disagreement between the flat cosmological model 
and the BBN predictions on the one side and the first Doppler peak posi- 
tion and the secondary peaks amplitudes in CMB power spectra obtained 
by BOOMERANG and MAXIMA on the other side. Negative values of few 
percent f or A n have b e en reported from BOOME RANG and MAXIMA data 



analysis flBattve et al I . l200ll ; lAvelino et al. [ l2000t) an d from BOOMERANG, 



DASI and COBE data analysis flAvelino et al. 1 . 120011 ). 



From the analysis of the fir st year WMAP data, a —0.06 < A a < 0.01 interval 



(IRocha et al. 1 . 12004 120031 ) was reported. Using the same data togheter with 



HST Hubble Key Project, but considering simultaneous time variations of 
both a and the electron mass m e , and using a different analysis method, 
has been obtained the interval —0.048 < A Q < 0.032 when only a varies 



2 



( Ilchikawa et al 



2006|). 



In the present work, the WMAP-3yr CMB data are being analyzed in order to 
find new limits on the a value at recombination, and the corresponding limits 
on the recombination redshiftf 1 "] 

In section 2 it is reviewed the standard recombination process and are de- 
scribed the changes involved by a different value of the fine structure constant 
at recombination. 

Section 3 presents the results of the analysis of WMAP-3yr data in the context 
of the recombination process affected by time varying fine structure constant. 

In section 4 are summarized the conclusions of the work. 



2 TIME VARYING a AND CMB 



It is well known the effect of changing the value of the fine structure con- 
stant on the energy levels of Hydrogen atom and, therefore, on the photon 
frequencies corresponding to the transitions between two such levels. As a 
consequence, a change in a value will change the dynamics of the Hydrogen 
recombination process in the Universe, which affects the CMB fluctuations 
observed today. As the CMB primary anisotropies probe the recombination 
epoch, it is natural to consider imposing limits on possible variation in a value 
at the recombination epoch relative to present value using CMB anisotropy 
data. 

In the following I briefly review the Hydrogen and Helium recombination pro- 
cess and the implications of non-standard value of a, neglecting the impact 
on Helium. 



Before recombination, the photon gas was coupled to electron-baryon fluid 
through Thoms on scattering on fr ee electrons, the cross section for this process 



being given by ([Weinberg I . Il995h 



1 e 4 , 

a T = oc a 

07T mi 



The CMB formation followed the recombination process, when the photons 
became free particles after electrons have been captured by ions. 



Throughout this paper, we consider the recombination redshift as being the red- 
shift corresponding to maximum of the visibility function. 
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For modeling the recom bination process I u sed the treatment implemented 



in the RECFAST code (ISeager et al. I . Il999l ). The equations describing the 



evolution of proton fraction x p , the singly ionized Helium fraction XHeii an d 
the matter temperature T M are: 



dz 

H(z)(l + z) [XeXpUHRH ~ Ph{1 ~ x v)^ hVH ' kTM l ( 2 ) 

dXHell _ Cue 

dz ~H(z)(l + z) 

x [x HeII x e n H R HeI - (3 HeI {f He - x HeII )e- h ^l kT % (3) 

dT M 



dz 



-{T M -T R ) + —^-. (4) 



3H(z)(l + z)m e l + f He + x e x nj (1 + z) 



In the above equations, H(z) is the Hubble expansion rate at redshift z, h 
is the Planck constant, k is the Boltzmann constant, c is the speed of light, 
or = fc 4 /(1207rc 3 /i 3 ) is the blackbody constant, nn is the Hydrogen number 
density, x p = n e /nn is the proton fraction, xueii = nneii/nn is the singly 
ionized Helium fraction and x e = n e jnn = x p + xueii is the electron fraction. 
The number fraction of Helium to Hydrogen is fn e — Y p /(4(1 — Y p )), where 
Y p = 0.24 is the primordial Helium mass fraction. The radiation temperature 
Tr(z) — TcmbO- + z) is identical to Tm at high redshift because of the cou- 
pling of photons and baryons through Thomson scattering. Rh is the case B 
recombination coefficient for H, and is given by the fit formula 

R H = F1(T 19 ^— mV 1 (5) 
l + ct d V ' 



where t = T A „/{W)K a = 4.309, b 



-0.6166, 



(IPequignot et al. I . ll99ll ) and F = 1.14 is the fudge factor (ISeager et al. I . ll999l ) 



0.6703, d = 0.5300 



introduced in order to reproduce the results of the multilevel calculation by 
speeding up recombination in the standard scenario. /3# is the photoionization 
coefficient 



R 



2irm P kTi 



M 



H 



exp(- 



B 



H2s - 



kT, 



(6) 



Al 
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and Ch is the Peebles reduction factor 

r = [1 + K H A H n H (l - x p )} 
H [1 + K H (A H + fa)n H (l - x p )\ ' U 



which accounts for the presence of non-thermal Ly-a resonance photons. In 
the above, Bh2s = 3.4eV is the binding energy in the 2s energy level, vh = 
(Bnis — Bn2s)/h is the Ly-a frequency, Ah is the rate of decay of the 2s 
excited state to the ground state via 2 photons, and Kh = c 3 / (StiuhH (z)). 

The quantities Rnei, @Heii C-ffe and i>Hei from Eq. (3) are the analogous for 
Helium of the quantities from Eq. (2) and their expressions may be found, for 
example, in Ichikawa et al. (2006). In this work, the small effect of changing 
a on Helium recombination process has been neglected. 

In order to take into account the changing value of a in the recombination 
process, I have modified the evolution equations for proton fraction and for 
matter temperature, considering their a dependence. The quantities which 
depends on a in these equations have bee n Taylor expanded up to first order 



in A a according to their scaling relations (jKaplinghat et al. I . Il999l ) 



Rh oc a 
B Hn oca 2 , 
Kh oc a -6 



2(1+0 



Ah oca 8 , (8) 
where £ = 0.7 was adopted. 



3 CMB CONSTRAINTS ON a VARIATION USING LATEST 
WMAP DATA 



In this work, in order to search for new CMB limits on a value at recombina- 



tion, I have analyzed th e WMAP-3yr CMB anisotropy data (IHinshaw et al 



20061 ; iPage et al. 1 . 120061 ) in the framework of the extended cosmological model 



which includes the variation of the fine structure constant at recombination 
with respect to its present value. The data anal ysis has been done using 



Mark ov Chain Monte Carlo (MCMC) techni ques (Isee e.g. MacKav. D.J.K 



jpoih implemented in the COSMOMC code flLewis fc Bridle I . l20Q2h 



For this purpose I modified the RECFAST code to compute the Hydrogen 
recombination in the hypothesis of different value of a at last scattering. The 



relative variation A r 



a )/a was added as an additional parameter; 
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the modified equations for the evolution of Hydrogen ioniza tion fraction and 



matt er temperature have been integrated with CAMB code (ILewis. A. et al. 



2000l ). used by COSMOMC to compute the theoretical CMB power spectra. 



Together with A a the following cosmological parameters have been varied: 
physical density in baryons Q^h 2 , physical density in DM Q^Mh 2 , Hubble 
constant H , reionization redshift z re , spectral index n s and amplitude A 2 R of 
primordial fluctuations; the cosmological constant Q\ and the optical depth 
to reionization r have been derived. To realistically const rain the Hubble ex - 
pansion rate taking into acco unt its degeneracy with A^, (IHannestad I . Il999i ) 



the HST Key Proje ct prior (IFreedman et al. I . l200ll ; ILewis fc Bridle I . 12002 



Rubino et al. 1 . 120031 ) has been used in addition to WMAP-3yr data. 



The modified version of COSMOMC software package has been run on 8 
Markov cha ins, using the "variance of chain means" /"mean of chain variances" 
R statistic (IBrooks fc Gelman I . Il998l ) as convergence criterion with the choise 
R- 1 < 0.03. 



The most likely values of cosmological parameters obtained from MCMC sim- 
ulations are given in Table 1. The mean values of the standard cosmologica l 



parameters are in the limits reported by WMAP team (ISpergel et al. 1 . 120061 ) . 



Table 1: ACDM with varying a model parameters and their 68% confidence 
intervals obtained from WMAP-3yr data. 



Parameter 


Mean 


ML 


won b h 2 


z - i0 -o.n 


2.18 


n DM h 2 


n 104+0-008 

U.1U^_q 00g 


0.108 


H 


68.43+jj 


68.74 




1 1 Q7+2.58 
LL -' i ' -2.56 


12.02 


n s 


Q6+ - 02 
u ' yu -0.02 


0.96 


W 10 A 2 R {k = 0.05) 


20.931^ 


21.39 


Age/GYr 




14.05 


A a 


-0 nil+ - 017 
u.uii_ 006 


-0.006 


Q A 


n 79 +0.06 
u - ' z -0.05 


0.73 




n 97+O.O6 
u ' z '-0.05 


0.27 


T 


0.090l°;°2 


0.096 



The marginalized distributions of the parameters obtained from simulations 
are presented in Figure 1. The observed differences between the curves corre- 
sponding to the parameters H , Q b h 2 , Q m , Q\ and the age of the Universe 
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Fig. 1. Marginalized distributions for the parameters of ACDM model with varying 
a obtained from the analysis of WMAP-3yr data and HST Key Project (continuous 
lines) compared with the distributions of the parameters in the standard ACDM 
model (dashed lines) from the same data. 



in the ACDM with varying a model compared to those of the standard ACDM 
model are due to the degeneracy of these parameters with a. 

Figure 2 presents the 2D marginalized distributions for pairs of parameters 
at 68% and 95% CL for ACDM with varying a and the ACDM models. The 
negative values for A Q in the ACDM with varying a model are correlated with 
smaller values of Hq and flbh 2 then in the ACDM model. On the other side, the 
distributions for Vt^Mh 2 are almost the same in both cases, therefore, in the 
case with varying a the parameter Q m = Q b +Q DM tends to go to higher values 
and to lower values then in the case with no a variations. Consequently, 
the distribution for the age of the Universe is widened, with mean value larger 
then in the case with constant a. 

Figure 3 shows the Cf T , Cf E and Cf E CMB power spectra for the best fit 
parameters of the model with varying a compared with the best fit CMB 
power spectra for the ACDM model and the WMAP-3yr experimental data. 
It is remarkable that the two best fit models are almost identical from the 
CMB power spectra point of view. 

The 95 % confidence interval for the additional parameter obtained from 
WMAP-3yr data is 
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Fig. 2. 2D - marginalized distributions for the parameters of ACDM model with 
varying a obtained from the analysis of WMAP-3yr data and HST Key Project 
(continuous lines) compared with the 2D - distributions of the parameters in the 
standard ACDM model from the same data (dashed lines). The contours are at 68% 
and 95% CL. 

(9) 

-0.039 < A Q < 0.010, 



which is a tighter CMB constraint on the variation of fine structure constant 
at rec ombination than previou s ly reported from the analysis of WMAP-lyr 



data (IRocha et al. 1 . 12004 120031 : llchikawa et al. 1 . 120061 ) . 



The interval for the recombination redshift corresponding to the limits on A c 



1 012 < z rec < 1 115, 
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Fig. 3. CMB power spectra C{ 



rr 



(a), Cf E (b) and Cf E 



c) for the best fit ACMD 



model with varying a (continous lines) and for the the best fit ACDM model wit h 



constant a (computed by the WMAP team (|lambda.gsfc.nasa.gov website I . |20o3 )) 



(dash-dotted lines) and the WMAP-3yr binned data and their errors. 



is tighter com pared to those correspond i ng; to limits on A„ fro m WMAP-lyr 
data analysis (IRocha et al. 1 . 120041 . 120031 ; llchikawa et al. 1 . 12006|). 



The recombination redshift corresponding to the best fit parameters of the 
model with varying a, z rec = 1 075, indicate a delayed recombination compared 
with the results from WMAP-lyr data. 

The relative variation of a in redshift unit corresponding to the best fit is 



a- l da/dz = -5.954 x lCT 6 
or, equivalent, in time unit, 



cT l da/dt = -4.65 x KT 13 ^ -1 , 



both in agreement with the pre dicted constraints on a from CMB (IHannestad 



19991 : iKaplinehat et al. I . Il999h . 
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4 CONCLUSIONS 



Based on the analysis of WMAP-3yr data with MCMC techniques, I obtained 
a tighter CMB constraint on the value of fine structure constant at recombi- 
nation epoch then was previously obtained from WMAP-lyr data. The results 
confirm the limits predicted to be achieved from CMB. 

The WMAP-3yr constraints on other cosmological parameters like Hq, Qbh 2 , 
fi m and fl\ are relaxed because of their degeneracy with the a variation, con- 
firming that a non-stand ard recombination weakens the constraints on other 



cosmological parameters (IBean et al. 1 . 120071 ) 



Using the constraint on a variation, I have also obtained the recombination 
redshift, finding a delay of the recombination epoch, similar to that of the 
reionization epoch and of the structur e formation computed from WMAP-3yr 



versus WMAP-lyr data (IPopa 1 . 120061 ). 
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